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ABSTRACT 

Jet-induced supernovae (SNe) have been suggested to occur in gamma-ray bursts (GRBs) and 
highly-energetic SNe (hypernovae) . I investigate hydrodynamical and nucleosynthetic properties of 
the jet-induced explosion of a population III 40Mq star with a two-dimensional special relativistic hy- 
drodynamical code. The abundance distribution after the explosion and the angular dependence of the 
yield are obtained for the models with high and low energy deposition rates -Edep = 120 x 10 51 ergs s _1 
and 1.5 x 10 51 ergs s _1 . The ejection of Fe-peak products and the fallback of unprocessed materials 
in the jet-induced SNe account for the abundance patterns of the extremely metal-poor (EMP) stars. 
It is also found that the peculiar abundance pattern of a Si-deficient metal-poor star HE 1424-0241 is 
reproduced by the angle-delimited yield for 9 — 30° — 35° of the model with Ed cp = 120 x 10 51 ergs s _1 . 
Furthermore, I compare the yield of the jet-induced explosion with that of the spherical explosion and 
confirm the ejection and fallback in the jet-induced explosion is almost equivalent to the "mixing- 
fallback" in spherical explosions. In contrast to the spherical models, however, the high-entropy 
environment is realized in the jet-induced explosion and thus [(Sc, Ti, V, Cr, Co, Zn)/Fe] are en- 
hanced. The enhancements of [Sc/Fe] and [Ti/Fe] improve agreements with the abundance patterns 
of the EMP stars. 

Subject headings: Galaxy: halo — gamma rays: bursts — nuclear reactions, nucleosynthesis, abun- 
dances — stars: abundances — stars: Population II — supernovae: general 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are phenomena emit- 
ting 7-ray for short periods followed by a power- 
law decaying afterglow. The origin had been cov- 
ered for a long while, but it has become clear 
that long GRBs are associated with supernovae 
(SNe). Three GRB-associated SNe h ave been ob- 
serve d so far: GRB 980425/SN 19 98bw (iGalama et all 
1998). GRB 030329/SN 2003dh (Stanek et al. 2003: 
Hiorth et all I2Q0J), and GRB 031203/SN 20031w 



(|Malesani et al.ll2~004[ ). They are all energ etic explosions 
of massive stars, called hypernovae (e.g.. [Nomoto et ahl 
120061 120071 for a review) . 

Although the explosion mechanism of GRBs and GRB- 
associated SNe is still uncovered, the following photomet- 
ric and spectroscopic observations indicate that they are 
asphcrical explosions with jet(s). (1) The light curve of 
the GRB afterglow has shown a polychromatic break in 
the power-law decay. The break is called a "jet break" 
explained by a deceleration of a relativistic jet and rel- 
ativis tic beaming of light (e.g., iFrail et aT1l2001t iPiran I 
120051 ). (2) The nebular spectra of SN 1998bw show nar- 
rower [O I] lines than [Fe II] lines (jPatat et al.l l2001h 
indicating that O locates in the inner and lower-velocity 
region than Fe. This is realized in a n aspherical explo- 
sion b ut not in a spherical explosion (|Maeda et alj|2002l 
I2006T ). 

The aspherical explosions are indirectly suggested 
from the abundance patterns of extremely metal-poor 
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(EMP) stars with [Fe/H] < -3. 3 Such EMP stars 
are likely to show the nucleosyn thesis yields of a sin- 
gle or a few core-collapse SN/SNe (|Audouze fe Silklll995l : 
iBeers fe Christliebl l2005f) . The abundance patterns are 
reproduced by the "mixing-fallback" models that as- 
sume the extensive mixing of the shoc ked material before 
the fallback in a spherical SN model (lUmeda fe Nomoto 
20021 120031 l2005Ulwamoto et alJl2005t iTominaga et al. 
2007b). In particular, the abundance patterns of EMP 
stars with high [Co/Fe] and [Zn/Fe] (~ 0.5) are explained 
by energetic explosion models because the high ratios re- 
quire explosive nucleosynthesis under high entro py (e.g., 
lUmeda fe Nomotol 120051 ITominaga et alj|2007b|) . Since 
the fallback doesnot take place in energetic spherical 
or quasi-spherical explosi ons (e.g., IWooslev fe Weaveil 
119951 : Hwamoto et al.ll2005l ). the mixing and fallback are 
interpreted as a consequence of the aspherical explo- 
sion. Indeed, the abundance patterns of the EMP 
stars are reproduced by the jet-induced explosion model 
(|Tominaga et al.ll2007al ). 

Since the jet-induced explosion contributes to many as- 
tronomical phenomena, it is important to make a quan- 
titative prediction on the nucleosynthesis outcome of the 
jet-induced explosions. In order to follow the jet prop- 
agation, fallback, and nucleosynthesis on the site and 
compare the jet-induced explosion models with the EMP 
stars, GRBs, and SNe, it is required to calculate multi- 
dimensional hydrodynamics and nucleosynthesis with the 
gravity and the special relativity. Although there are 
many studies on the GRB jet propagation of the stellar 

3 Here [A/B] = log 10 (iV A /iVB) - log lo (JV A /Ar B ) , where the 
subscript refers to the solar value and TVa and A?b are the abun- 
dances of elements A and B, respectively. 

4 In order to follow the core collapse, jet formation, and nucle- 
osynthesis on the site, it is also required to involve the neutrino, the 



mantle using the special relatiyistic hydrodynamics (e.g . , 
lAloy et al.ll2000HZhang et al.ll2004t iMizuta et alll2006h . 
they did not include the gravity or calculate nucleosyn- 
thesis. 

To investigate the yields of the jet-induced SNe, it 
is crucial to include the gravity because the fallback 
plays an important role on the nucleosynthesis yields. 
The studies by the use of Newtonian calculations have 
concluded that the energy depositio n rate (-Edep) sen- 
sitive l y affects SN nucleo synthesis (jMaeda fc Nomotol 
120031 : iNagataki et al.ll2006Q . This resul t has been con- 
firmed for the special rel atiyistic cases ([T bminag a et al.l 
l2007al) . In particular, iTominaga et all (|2007af) have 
shown that the jet-induced explosions with various E^op 
can explain both the variations of GRB-associated SNe 
and the EMP, carbon-enhanced EMP (CEMP, [C/Fe] > 
1), and hyper metal-poor (HMP, [Fe/H] < -5) stars. 

The previous studies have proved mostly the angle- 
integrated yields and shown that the abundance pat- 
terns of the EMP stars are reproduced by the angle- 
integrated yield. However, the abundance distribution of 
the jet-induced explosio n depends on the direction (e.g., 
iMaeda fc Nomotol I2003D . Thus, the abundance patterns 
of the next-generation stars might depend on the direc- 
tion. I calculate aspherical stellar explosions induced by 
highly relativistic jets and obtain hydrodynamical and 
nucleosynthetic structures of such jet-induced explosion 
models. In particular, I investigate the angular depen- 
dence of the yield to compare the yields with the abun- 
dance patterns of the metal-poor stars. Furthermore, I 
compare the jet-induced explosion with the spherical SN 
model applied the mixing-fallback model and connect 
properties of the jet-induced explosion to the mixing- 
fallback model. 

In § |U the applied models are described. In § 02 I 
present the hydrodynamical and nucleosynthetic struc- 
tures of the jet-induced explosion model, investigate the 
angular dependence of the yields, and compared the jet- 
induced explosion model with the spherical SN model. In 
§ |H the conclusion and discussion are presented. In Ap- 
pendixes, the hydrodynamics and nucleosynthesis code 
is described and tested. 

2. MODELS 

I investigate a jet-induced SN exp l osion of a Pop 
III40M Q star (jUmeda fc Nomotoll2005t ITominaga etTal] 
2007b) by means of a two-dimensional relativistic Eu- 
lerian hydrodynamic and nucleosynthesis calculation 
with the gravity (Appendix |A"]) . The nucleosyn- 
thesis calculation is performe d as a post-processing 
(jHix fc Thielemann! 1 19961 Il999h with the reaction net- 
work incl uding 280 isotopes up to 79 Br (see Ta- 
ble 1 in lUmeda k. Nomotol l2005h . The thermody- 
namic histories are traced by make r particles represent- 
ing Lagrangian mass ele ments (e.g.. lHachisu et al"1ll990l : 
IMaeda fc Nomotoll2003l see also Appendix[Bj . A compu- 
tational domain initially ranges up to the stellar surface 
where i? s tar = 2 x 10 7 km and is captured by 200 logarith- 
mical grids in the r-direction and 100 uniform grids in 
the 0-direction. A circumstellar matter (CSM) extends 
from the stellar surface with the slope p oc r~ 2 . 




Fig. 1. — Schematic picture of the jet-induced explosion. 



The explosion mechanism of GRB-associate d SNe is 
still under debate (e.g., a neut r ino an nihilation, IWooslevi 
119931: iMacFadven fc Wooslevl 119991 : and a magneto- 
rotation. iBland ford & Znaiekl Il977t iBrown et all [20001 : 
iMizuno et al.ll2004f ). Thus. I do not consider how the jet 
is launched, but the jet is dealt parametrically with the 
following five parameters (Fig. [I]): energy deposition rate 
(-Edep), total deposited energy (E'dcp), initial half angle 
of the jets (#j e t), initial Lorentz factor (Fj Ct ), and the ra- 
tio of thermal to total deposited energies (/th)- The jet 
is injected from the inner boundary at an enclosed mass 
Mq corresponding to a radius Rq. The density, pressure, 
and velocities of the jet are described with the five pa- 
rameters (Appendix [Cj and put as a boundary condition 
at the inner boundary. The jet is assumed to consist of 
the accreted matter and to expand adiabatically below 
the inner boundary. After the jet is injected into the 
computational domain, the thermodynamic histories are 
traced by the marker particles. 

In this paper, I show three models; (A) a model with 
E dcpM = ^dc P /(10 51 ergs s" 1 ) = 120 and M = 1.4M© 
(R = 900 km), (B) a model with E dcpM = 1-5 
and M = 1.4M (R = 900 km), and (C) a model 
with E dcp . 51 = 120 and M = 2.3M (R = 2700 
km). The other parameters are same for each model; 
E dcp = 1-5 x 10 52 ergs, 5 jot = 15°, r jct = 100 and 
/th = 10~ 3 . The mass of jets is M jet - 8 x 10- 5 M Q . 



general relativity, and possibly the magnetic field (e.g.. |Pruet et aLl 
2005; Fro hlich et aLl l2006: Janka et al. 2007|). 



5 IFrail et al.l (120011') suggested the 7-ray energies (E 7 ) of GRBs 
are clustered at E 7 ~ 5 X 10 50 ergs. Although E 7 is 30 times 
smaller than -Edep i there are large uncertainties on the relation be- 
tween E~f and -Edep- F° r example, the energy possessed by the rela- 
tivistic matter depends on the interactions between the relativistic 
jet and the stellar mantle and is reduced compared to the injected 
energy (Zhang et al. 2003, 2004). Moreover, the radiative conver- 
sion efficiency of a kinetic energy depends on an unknown 7-ray 
emission mechanism. The efficie ncy is estimated to be ~ 2 — 40% 
for collisions of internal shocks ( Kobayashi et al. 1997) or > 50% 
for nonthermalized photospheric emissions (Ioka et al. 2007). -Edep 
adopted in this paper might be appropriate for the former emission 
mechanism or excessive for the latter emission mechanism. The de- 
pendence of the jet-induced explosions on -Edep will be studied in 
future. 



Name 



TABLE 1 

Jet-induced explosion models. 



M 
[M e ] 



Ro 
[km] [10 51 ergs s" 1 ] [10 51 ergs] 



-^dcp 



^clcp 



Pjct 

[degrees] 



M rom 
[M Q ] 



A 


1.4 


900 


120 


15 


15 


100 


10~ 3 


9.1 


B 


1.4 


900 


1.5 


15 


15 


100 


10" 3 


16.9 


C 


2.3 


2700 


120 


15 


15 


100 


10" 3 


8.1 




Fig. 2. — Density structures at t = 10 5 s along the jet axis (solid 
line) and the equatorial plane (dashed line) of model A (thick line) 
and B (thin line). 

The model parameters and the central remnant mass 
(M rcm ) are summarized in Table [TJ Models A and B 
are used in iTominaga et al.l (|2007al ) and th ey reproduce 
the a bundance patterns of the EMP stars 6 (ICavrel et al 
200l and HE 13 27-2326 (e.g., iFrebel et al] T2005. 200 
Collet et al.ll2006l ). respectively. 



3. RESULTS 

The hydrodynamical calculations are followed until the 
homologously expanding structure is reached (v ex r). 
Then, the ejected mass elements are identified from 
whether their radial velocities exceed the escape veloci- 
ties at their positions. The density structures of models 
A and B at t = 10 5 s are shown in Figure^ The den- 
sity along the jet axis is higher than the density along the 
equatorial plane because the matter is more easily ejected 
along the jet axis than along the equatorial plane. And 
the SN ejecta of model A is denser and more compact 
than that of model B because the ejected mass of model 
A is larger than that of model B. 

3.1. Fallback 

Figures [3^ and [3Jd show "accreted" regions for mod- 
els A and B, where the accreted mass elements initially 
located in the progenitor. The O layer is separated into 
the two layers: (1) the O+Mg layer with X( 24 Mg) > 0.01 
and (2) the O+C layer with X( 12 C) > 0.1. The inner 
matter is ejected along the jet-axis but not along the 

6 An averaged abundance pattern of four EMP stars, CS 22189- 
009, CS 22172-002, CS 22885-096, and CD -38 245, is adopted. 
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Fig. 3. — Initial locations of the mass elements which are finally 
accreted (black), for (a) model A and (b) model B. The background 
circles represent the boundaries between the layers in the progeni- 
tor star; the H, He, O+C, O+Mg, and Si layers from the outside. 

equatorial plane. On the other hand, the outer matter is 
ejected even along the equatorial plane, since the lateral 
expansion of the shock terminates the infall as the shock 
reaches the equatorial plane. 

The accreted mass is larger for lower -Edcp- This stems 
from the balance between the ram pressures of the i nject- 
ing jet (Pjet) and the infalling matter (Pf a ii) (e.g., iFrverl 



13; tv 



[J99a [Fryer - fc Meszarost2003l:lMaeda fc Tominagal2008ft . 
While Pj t depends on i^dop, -Pfaii depends on the presu- 
pernova structure. If the matter falls freely, Pf a n is pro- 
portional to p r r 3 / 2 M r for the matter having located at r 



and M r in the progenitor where the presupernova density 
is p r . The critical energy deposition rate (-Edep.cri) giving 

-Fjet (^dep.cr i) = Pfai\(M r ) is lower for the outer layer (i.e., 
large r M r , iFrver fe Meszarosl [200l iMaeda fe Tominagal 
2008). Therefore, the jet injection with lower Ed ep is en- 
abled at a later time when the central remnant becomes 
more massive. Additionally, the lateral expansion of the 
jet is more efficiently suppressed for lower -Edcp- As a 
result, the accreted region and M rem are larger for lower 

Edcp- 

A model with lower -Edcp has larger M rcm , higher 
[C/O], [C/Mg], and [C/Fe], and smaller M(Fe) becaus e 
of the larger amount of fallback (|Tominaga et al.ll2007af ). 
The larger amount of fallback decreases the mass of the 
inner core (e.g., Fe, Mg, and O) relative to the mass of 
the outer layer (e.g., C, Figs.[3kb). Since O and Mg are 
synthesized in the inner layers than C, [C/O] and [C/Mg] 
are larger for the larger infall of the layer. Also, the 
fallback of the O layer decreases M(Fe) because Fe is 
mainly synthesized explosively in the Si and O+Mg lay- 
ers. Therefore, the variation of £dcp in the jet-induced 
explosions predicts that the variations of [C/O], [C/Mg], 
and [C/Fe] are corresponding to the variation of M (Fe). 
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3.2. Abundance distribution 

Figures 0^ and 0Jd show the abundance distributions 
and density structures at t = 10 5 s for models A and B. I 
classify the mass elements by their abundances as follows: 
(1) Fe with X( 56 Ni) > 0.04, (2) Si with A( 28 Si) > 0.08, 

(3) O+Mg with A( 16 0) > 0.6 and A( 24 Mg) > 0.01, 

(4) O+C with X( lh O) > 0.6 and A( 12 C) > 0.1, (5) 
He with A( 4 He) > 0.7, and (6) H with X( l R) > 0.3. 
If a mass element satisfies two or more conditions, the 
mass element is classified into the class with the smallest 
number. 

The abundance distribution and thus the composition 
of the ejecta depend on the direction. In model A, the 
O+Mg, O+C, He, and H mass elements locate in the all 
direction. On the other hand, most of the Fe and Si mass 
elements locate at 9 < 10° and stratify in this order from 
the jet axis and a part of them locate at 15° < 9 < 35°. 
Interestingly the Fe mass elements surround the Si mass 
elements at 15° < 9 < 35°. In model B, most of the 
O+C and He mass elements locate at 9 < 3° because 
they are ejected only along the jet axis, while the Fe 
mass elements injected as a jet expand laterally up to 
9 ~ 50° and the H mass elements are distributed in the 
all directions. The lateral spread of the Fe mass elements 
in models A and B are led by the collision with the stellar 
mantle and the internal pressure of the jet. 

3.2.1. Angular dependence of the yields 

I investigate the angle-delimited yields integrated over 
each 10° for models A and B, although the integra- 
tion range might be too wide to be taken in a sin- 
gle next-generation star. Fe is mainly distributed at 
9 < 40° for model A (Fig. UK) and 9 < 50° for model 
B (Fig.[3J)). Figures [5^ andEb show the abundance pat- 
terns of the angle-delimited yields for 9 < 40° of model 
A and 9 < 50° of model B, respectively. The yields are 
compared with the abundance patterns of the EMP stars 
(|Cavrelet all 120041. the CEMP stars (HE 1300+0157, 
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Fig. 4. — Density structure (back ground gray scale) and the 
positions of the mass elements at t = 10 5 s for (a) model A and 
(b) model B. Color of the marks represents the abundance of the 
mass element (H: yellow, He: cyan, O+C: green, O+Mg: blue, Si: 
magenta, and Fe: red). Size of the marks represents the origin of 
the mass element (the jet: dots, and the shocked stellar mantle: 
filled circles). 



IFrebel et alJl2007h. an d the HMP stars (HE 0107-5240 
e.g.. IChristlieb et all l2002t iBessell fe Christliebl 120051: 
i Collet et all 120061 and HE 13 27-2326. e.g.. IFr ebel et all 
2005, 2005" lCollet et all 12001) . 

FigureEK shows the angle-delimited yields of model A. 
The abundance patterns of the angle-delimited yields are 
determined by which mass elements are included into the 
integration. Because of the stratified abundance distri- 
bution, the yields for 0° < 6 < 10° and 10° < 9 < 20° 
show low [C/Fe] and [O/Fe] and high [C/Fe] and [O/Fe], 
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Fig. 5. — Comparisons between the abundance patterns of 
the angle-integrated yiel d (black), the EM P, CEMP, and HMP 
stars (red filled marks, ICavrel et alj I2004I; IFrebel et al.l I 2007I ; 
Christlieb et al. 2002; Bcsscl l fe Christlieb|2 005; Colle t "eTaIj20od ; 
Frebel et alj |2005, 200l|), and the angle-delimited yields (colored 
lines) for (a) model A and (b) model B. The color of the lines 
represent the yields integrated over 0° < 9 < 10° (green line), 
10° < 9 < 20° (blue line), 20° < 6 < 30° (magenta line), 
30° < 9 < 40° (cyan line), and 40° < 9 < 50° (gray line), (c) 
Comparison between the abund ance pattern of HE 1424-0241 (red 
filled circles, Cohen et al. 2007) and the angle-delimited yields of 
model A for 30° < 9 < 40° (solid line) and 30° < 9 < 35° (dashed 
line) . 

respectively. Intriguingly, the region with 30° < 9 < 40° 
includes the Fe mass elements in the outer layer and the 
O+Mg and O+C mass elements in the inner layer but not 
the Si mass elements. Thus, the yield for 30° < 9 < 40° 
shows a Si-deficient abundance pattern. Furthermore, 
the Fe and Si mass elements at 15° < 9 < 35° initially 
located along the jet axis (9 < 10°), and thus the high- 
entropy environment is realized in these mass elements 
and [Sc/Fe], [Ti/Fe], [Co/Fe], and [Zn/Fe] are enhanced 
in the yields for 10° < 9 < 40°. On the other hand, 
[Sc/Fe] is low in the yield for 0° < 9 < 10° because the 
matter initially located at 10° < 9. 

Figure Efc> shows the angle-delimited yields of model 
B. The all angle-delimited yields of model B show high 
[C/Fe]. This is the same feature as the angle-integrated 
yield. Since most of the heavy elements locate at 9 < 10°, 
the abundance pattern of the yield for 0° < 9 < 10° is 
similar to that of the angle-integrated yield, except for N 
that is mostly contained in the H mass elements. On the 
other hand, the yields for 10° < 9 < 50° consist of the 
Fe mass elements injected as a jet and the O+C mass 
elements. In the yields for 10° < 9 < 50°, [C/Mg] is 



almost similar and the differences of [C/Fe] and [Mg/Fe] 
mainly stem from the different amount of Fe. 

3.2.2. Abundance patterns of the metal-poor stars 

A very peculiar, Si-d eficient, metal-poor star HE 1424- 
0241 was observed (jCohen et al.1 [2007). Its abun- 
dance pattern with high [Mg/Si] (~ 1.4) and nor- 
mal [Mg/Fe] (~ 0.4) is difficult to be repro- 
duced by previous SN models. This is because 
log{[X(Mg)/X(Si)]/[X(Mg)/X(Si)] } < 1.6 is real- 
ized in the O+Mg layer at the presupernova stag e 
(e.g.. lWooslev fe Weaverll 19951: lUmeda fe Nornoto| [2005). 
Thus, in order to reproduce the abundance pattern 
of HE 1424-0241, the SN yield is required to in- 
clude explosively-synthesized Fc but not explosively- 
synthesized Si. 

The angle-delimited yield may possibly explain the 
high [Mg/Si] and normal [Mg/Fe] (Fig. Gib). Figure Et 
shows that the yields integrated over 30° < 9 < 40° and 
30° < 9 < 35° of model A reproduce the abundance 
pattern of HE 1424-0241. The yield consist of Mg in 
the inner region and Fe in the outer region (Fig. Hk). 
Although there are some elements to be improved, the 
elusive feature of HE 1424-0241 could be explained by 
taking into account the angular dependence of the yield. 
The high [Mg/Si] and normal [Mg/Fe] can be realized 
with an appropriate integration range if the Fe mass el- 
ements penetrate the stellar mantle (i.e., the duration of 
the jet injection is long) and if the O+Mg mass elements 

are ejected in all directions (i.e., -B^ep i s high). 

If the yield depends on the direction, the abundance 
patterns of the angle-delimited yields have a large scatter 
(Fig. [5k.). Although the diversities of [C/Fe] and [O/Fe] 
in the angle-delimited yields are similar to those in the 
observations of the EMP and CEMP stars, the observa- 
tions of the EMP stars provide abundanc e ratios with 
comp aratively small scatters (e.g., [Sc/Fe] ■ fCavrel et al.1 
2004) . This implies that SNe contributing to most of the 
EMP stars experienced strong mixing of the ejecta, e.g., 
due to an interaction with interstellar medium (ISM), 
so that the yield no longer depends on the direction. 
Actually, it was suggested that the interaction between 
the SN ejecta and ISM induces further mixing of the 
abundances if ISM has an inh omogeneous structure (e.g., 
iNakasato fe Shigevamal 120001 ) . According to stru cture 
formation calculations (e.g., lYoshida et al.l 120071 ). the 
Pop III stars are surrounded by the inhomogeneous ISM. 
Therefore, the mixing of the SN ejecta is likely to take 
place. On the other hand, the abundance pattern of 
HE 1424-0241 implies weak mixing of the SN ejecta. 
Such a different star formation process might cause the 
peculiar abundance pattern of HE 1424-0241. In order 
to conclude the origin of HE 1424-0241, it is necessary 
to calculate three-dimensional evolution of a supernova 
remnant in the inhomogeneous ISM. 

3.3. Comparison with the spherical supernova model 

The calculations of the jet-induced explosions show 
that the ejection of the inner matter is compatible with 
the fallback of the outer matter (Figs. [3kb). This 
is consistent with the two-dimensional illustration of 
the m ixing-fallback model (Fig. 12b in iTominaga et al.l 
I2007H) . 
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Fig. 6. — (a) Density structure (back ground gray scale) and the positions of the mass elements at t = 10 s for model C. The color and 
size of the marks represent the same as Figs. |4^,b. (b) Initial locations of the mass elements which are finally accreted (black) and ejected 
(colored points) for model C. The color of the mass el ements rep resents the same as Figs.[4^,b. (c) Comparison of the abundance patterns 
of model C (black line), the EMP stars (red circles, Cayrcl ct al. 2004), and the mixing-fallback models with (1) M cu t(ini) = 2.3Mq, 
M mix (out) = 12.2M© and / = 0.41 (green line), (2) M cut (ini) = 2.3M Q , M mix (out) = 10.3M Q and / = 0.27 (blue line), and (3) 
Mcut(ini) = 2.3M Q , M mix (out) = 10.8M© and / = 0.19 (magenta line). 



In this subsection, I clarify the relation between the 
mixing-fallback model and the jet-induced explosion 
model by comparing the yields. The mixing- fallback 
model has three parameters; initial mass cut [M cu t(ini)], 
outer boundary of the mixing region [M m ; x (out)], and a 
fraction of matter ejected from the mixing region (/). 
The remnant mass is written as 

M rom = M cut (ini) + (1 - /)[M mix (out) - M cut (ini)]. (1) 

The three parameters would relate to the hydrodynam- 
ical properties of the jet-induced explosion models, e.g., 
the inner boundary (Mo), the outer edge of the accreted 
region (M acc , out ), and the width between the edge of the 
accreted region and the jet axis. 

I apply model C for the comparison. The inner bound- 
ary, the outer edge of the accreted region, and the 
central remnant mass of model C are Mq = 2.3M Q , 



M, 



acc,out 



12.2Af Q , and M r , 



8.1M m . The abun- 



dance distribution of model C at t = 10 s is shown in 
Figure [5h. 
Figure[6j3 shows the initial positions of the ejected mass 



elements in model C. Color of the marks represents the 
composition of the mass elements. Explosive nucleosyn- 
thesis takes place in the mass elements between the edge 
of the accreted region and the jet axis. Therefore, the 
amounts of the explosive nucleosynthesis products are 
reduced relative to the spherical explosions. The shock 
collision at the equatorial plane realizes the mass ejection 
from the deep O+Mg layer along the equatorial plane. 
Although the mass ejection could be due to the reflective 
boundary condition on the equatorial plane, this does not 
affect the following discussion because the amount of the 
mass ejection is small. 

The angle-integrated yield of model C is compared with 
the yields of the spherical SN model with M ms = 40M© 
and £dcp = 3 x 10 52 ergs. For the spherical SN model, 
the explosion energy is deposited instantaneously as a 
thermal bomb. Here, I set M cut (ini) to be the same as 
M [i.e., M cut (ini) = 2.3M Q ] and apply three sets of 
parameters, M m j x (out) and /, as follows: 

(1) M mix (out) = 12.2M and / = 0.41. M mix (out) 
corresponds to M accout and / is the ratio of the ejected 



mass and the total mass between Mq and M acC)OU t. 

(2) M mix (out) = 10.3M© and / = 0.27. / is the frac- 
tion of the solid angle of the region where the ejected 
mass elements located at r ~ 10 3 km (the Si-burning re- 
gion, hereafter the fraction is written as /sO- Af m j x (out) 
is set to yield the same M rcm as model C. 

(3) M mix (out) = 10.8M© and / = 0.19. / and 
Af m i x (out) are set to yield the most resemblant abun- 
dance pattern to model C. The resultant M rem (= 
9.2M©) is slightly larger than M rem of model C. 

Figure [6J; shows that the angle-integrated abundance 
pattern of model C is roughly reproduced by the mixing- 
fallback model with the parameter sets of (2) and (3). 
Thus, I conclude that the jet-induced explosion is ap- 
proximated by the mixing-fallback model reasonably. 
Mq, /si, and M rcm in the jet-induced explosion model 
are represented by Af cut (ini), /, and M rcm in the mixing- 
fallback model, respectively. 

There are some elements showing differences, Sc, Ti, 
V, Cr, Co, and Zn. The enhancements of [Sc/Fe] and 
[Ti/Fe] improve agreements with the observations. The 
differences stem from the high-entropy explosion due 
to the concentration o f the energy injection in the jet- 
induced explosion (e.g-. lMaeda fe Nomotoll2003T) . In par- 
ticular, Sc is the most sensitive element to the entropy 
and [Sc/Fe] is more enhanced for the deeper Mq because 
of the weaker lateral expansion (Fig. GJi). Such thermo- 
dynamical features of the jet-induced explosion model 
cannot be reproduced by the mixing-fallback model ex- 
actly but a "low-density" modification might mimic the 
high-entropy environm ent fe.g.. iUmeda fe Nomotoll2005t 
iTominagaet alj|2007b| ). 

4. CONCLUSIONS AND DISCUSSION 

I perform two-dimensional hydrodynamical and nucle- 
osynthesis calculations of the jet-induced explosions of a 
Pop III 40M© star. Here I test three jet-induced explo- 
sion models A, B, and C as summarized in Tabled] and 
conclude as follows. 

(1) Fallback: The dynamics and the abundance distri- 
butions depend sensitively on the energy deposition rate 
-Edcp- The explosion with lower -Edep leads to a larger 
amount of fallback, and consequently smaller M(Fe) and 
higher [C/O], [C/Mg], and [C/Fe]. Such dependences of 
[C/Fe] and Af(Fe) on E dcp predict that higher [C/Fc] 
tends to be realized for lower [Fe/H] . This is consisten t 
with the observations (e.g., iBeers fe Christliea 120051) . 
Note, however, the formation of star with low [C/Fe] and 
[Fe/H] is possible because [Fe/H] depends on the swept- 
up H mass, i. e., the int eraction between the SN ejecta 
and ISM (e.g., ICioffi et al.lfl98l . 

(2) Abundance distribution and angular depen- 
dence: I present the aspherical abundance distributions 
and investigate the angular dependence of the yield. Af- 
ter the explosion, the elements ejected even along the 
equatorial plane locate in the all directions, while the el- 
ements ejected only along the jet axis expand laterally 
but tend to be distributed along the jet axis. The abun- 
dance distributions in the jet-induced SN ejecta could be 
examined by spatially-resolved observatio ns of supernov a 
remnants fe.g. JHwang et aL 2004; Kcoh ane et al.ll2007f ). 

The angle-delimited yield could reproduce the ex- 
tremely peculiar abundance pattern of HE 1424-0241. 
However, the angle-delimited yields of model A have a 



large scatter that may be inconsistent with the relatively 
small scatter in the abundance ratios of the EMP stars. 
This implies that the angular dependence of the yield 
in most SNe is diluted by the strong mixing of the SN 
ejecta. On the other hand, the angle-delimited yields of 
model B show high [C/Fe] that is the same feature as the 
angle-integrated yield. 

The angle-delimited yield strongly depends on which 
mass elements are included into the integration. This 
would be determined by the abundance mixing in the 
SN ejecta and by the region where the next-generation 
star takes in the metal-enriched gas. To investigate this 
issue further, it is required to calculate three-dimensional 
evolution of the supernova remnant in the ISM. 

(3) Comparison with the spherical explosion: 
The angle-integrated yield of the jet-induced explosion 
is well reproduced by a spherical SN model applied the 
mixing-fallback model. This confirms that the mixing- 
fallback model approximates the jet-induced explosion 
well and that the mixing and fallback in hypernovae as- 
sumed in the mixing- fallback model are actually achieved 
in aspherical explosions. The abundance ratios between 
elements synthesized in different regions (e.g., C, O, Mg, 
and Fe) depend on the hydrodynamical structure of the 
explosion, e.g., the fallback. Thus, such macroscopic 
properties of the jet-induced explosion are represented 
by the mixing-fallback model. In particular, Mo, /si, 
and M rem in the jet-induced explosion model are repre- 
sented by Mcut(mi), /, and M rom in the mixing- fallback 
model, respectively. 

On the other hand, the ratios between the explosively- 
synthesized elements depend on the thermodynamical 
properties of the explosion. In particular, [Sc/Fe], 
[Ti/Fe], [V/Fe], [Cr/Fe], [Co/Fe] and [Zn/Fe] are en- 
hanced by the high-entropy environment in the jet- 
induced explosion, thus showing differences from the 
mixing-fallback model. The enhancement of [Sc/Fe] 
and [Ti/Fe] improve the agreement with the observa- 
tions. Note, the enhancement of [Ti/Fe] relative to 
the spherical SN model in this paper is larger than in 
iMaeda fe Nomotol (|2003l ). This might be because the 
relativistic jet suppresses the lateral expansion and en- 
hances the energy concentration. 

The enhancement of [Sc/Fe] is also suggested to be 
obtained by nucleosynthesis in the p-rich ejec t a ex- 
pelled from the innermo s t region dPruet et al.l 120041 
120051 iFrohlich et all 120061 : llwamoto et alJl2006f ). How- 
ever, it seems difficult to enhance [Zn/Fe] (~ 0.5) and 
[Sc/Fe] (~ 0) on ly with nucleosynthesis in the p-rich 
ejecta (see § 5.2 in Tominaga et al. 2007b). On the other 
hand, a jet-induced explosion enhances [Sc/Fe], [Ti/Fe], 
and [Zn/Fe] simultaneously and realizes the mixing and 
fallback in an energetic explosion. Therefore, I propose 
that the jet-induced explosion is a presumable origin of 
the enhancement of [Sc/Fc]. The origin would be con- 
cluded by future quantitative investigations including 
various elements. 

Although the models with same #j e t are shown in this 
paper, the asphericity is likely to be different for each 
SN. For example, the observations of supernova rem- 
na nts show various aspherical structures (e .g., Cassiopeia 
A. iHwang et alJ [20 04: Fe sen et all 120061 : and W49B, 
iMiceli et al.l l2006t iKeohane et al.l 12007). Furthermore, 
recent observations of the nebular spectra of SNe sug- 



gested that all core-colla pse SNe are aspher ical explo- 
sions (jMaeda et alJ l2008t iModjaz et al J 12008). The de- 
gree of asphericity in SNe in the early universe will be 
revealed by detailed comparisons of abundance ratios be- 
tween the aspherical SN models with different 0j Ct and 
the metal-poor stars which represent the hydrodynami- 
cal and/or thermodynamical properties. 
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APPENDIX 
SPECIAL RELATIVISTIC HYDRODYNAMIC CODE 

A two-dimension al special relativistic Eulerian hydrodynamic code is developed with Marquin a's flux formula 
(iDonat et al.|[l998l ) and with a conversion method from the observer frame to the proper frame dMartf fc Midler 
119961 ). The code applies third order Runge-Kutta method in tim e of Shu fc Oshe r (1988, see also lAlov et al.lll999h 
and second order piecewise hyperbolic method (PHM) in space of lMarquinal (|1994h . 

The equations of special relativistic hydrodyn amics are described in terms of a four-velocity vector held and an 
energy momentum tensor (Mar ti fc Mullerlll994f ). Physical quantities in a rest frame (relativistic rest-mass density: 
-D, the i-th components of momentum densities: Si, and energy density: r) and a comoving frame are related as 
follows: 

Si=phT 2 Vi, 
t = phY 2 — p — pT, 
where the light velocity is set to c = 1, p is the proper rest-mass density, Vi are the i-th components of velocities, 
T = 1/yl — J2i v f 1S tne Lorentz factor of the fluid element with respect to the rest frame, p is the proper pressure, 
and h = 1 + e pr ,j n /p + pj p is the specific enthalpy (here, e pr] i n is the proper internal energy per unit volume). 
The basic equations in the special relativistic hydrodynamics are written in spherical polar coordinates (r, 9) as 



dD 

~dt 
dS r 
~df 

dS_ e , d(r 2 S g v r ) 

dt 



d(r 2 Dv r ) ld(sm9Dv e ) 



0(r 3 /3) 
d(S r v r + p) 
dr 



r d(—cos9) 
1 d(sm0S r vg) 
r d(-cosO) 
1 d(S e v e 



= 0, 



= -2- 



Of 



Sg vg 



■P) 



1 COS 9 SgV r 



dr 



<9(r 3 /3) r 89 

d{r 2 (S r - Dv r )} 1 d{sm9{S e - Dvg)} 



9rP, 



+ 9ep, 



cost 



: (g r v r + ggvg)p, 



(Al) 



(A2) 



(A3) 



(A4) 



<9(r 3 /3) r 0(- 

where gi (i — r, 9) are the i-th gravitational acceleration components. The equations are an equation of continuity, 
momentum conservation equations, an energy conservation equation. The gravitational potential includes the contri- 
butions of the self-gravity and the central remnant. I test two methods to include the self-gravity as follows: (1) the 
Poisson equation is approximated with an integration of spherical harmonics (Hachisu 1986) and (2) the sparse banded 
matrix r esulting from d ifferencing the Poisson equation is solved with a bi-conjugate gradient stabilized (BiCGSTAB) 
method (|Barrett et al.|[l994l ) 7 with modified incomplete LU (MILU) factorization. Since the results are consistent, I 
apply the former method in the calculations. 

One- dimensional Riemann problems 

The code is tested with one-dimensional shock tube problems (e.g. jMartf fc Mullerlll9941) . I set the computational 
region < x < 1 where c = 1 that is divided into two regions, i.e., left and right regions, at x = 0.5. The density (p), 
pressure (p) , and velocity (v) of left and right regions are represented b y L and R subscripts, respectively. The time 
variation is derived by the analytical solution of the Riemann problem (Ma rti fc Muller| [i994). I test three problems 
with I0 3 uniform meshes (Ax = 1/10 3 ) and a constant adiabatic index. The initial states are as follows: 

Problem (a) 7 = 5/3 

(pl,Pl,*'l) = (M0 3 ,0) at 0.0 < x < 0.5, 
{PR,PR, vr) = (1, 0.01, 0) at 0.5 < x < 1.0. 
Problem (b) 7 = 4/3 

(p L ,p L , T L ) = (1, 1, 10 3 ) at 0.0 < x < 0.5, 
(PR,PR, r R ) = (1, 10 6 , 1) at 0.5 < x < 1.0. 



Subroutines are available at Netlib (http://www.netlib.org/ ). 



Problem (c) 7 = 4/3 

(PUPL, vl) = (1, 10, -0.9) at 0.0 < x < 0.5, 
(PR, PR, vr) = (10, 100, 0.9) at 0.5 < x < 1.0. 

The analytical solutions are reproduced by the calculations (Figs. [ATk bc). The agreements confirm that the code 
correctly solves left- and right-oriented rarefaction and shock waves. 
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Fig. Al. — Comparison between the analytical solutions of the Riemann problems (black solid line) and the results of the special 
relativistic hydrodynamics calculations (density: red marks, pressure: green marks, velocity: magenta marks, Lorentz factor: cyan marks). 



Two-dimensional shock tube problem 

The two-dimensi onal shock tube problem was proposed by iDel Zanna fc Bucciantini (12002T) and repeated by subse- 
quent studies (e.g.. lLucas-Serrano et al.ir2004HZhang fc MacFadvenll2006t iMizuta et al.ll2006l ). I test the problem with 
10 3 x 10 3 square uniform meshes. In this problem, the two-dimensional region (0 < x < 1,0 < y < 1) are divided into 
the following four regions: 

(p,p, v x , v y ) = (0.50, 1.00, 0.00, 0.00) at 0.0 < x < 0.5, 0.0 < y < 0.5, 
(p,p, v x ,Vy) = (0.10, 1.00, 0.99, 0.00) at 0.0 < x < 0.5, 0.5 < y < 1.0, 
(p,p, v x ,v y ) = (0.10, 0.01, 0.00, 0.00) at 0.5 < x < 1.0, 0.0 < y < 0.5, 
{p,p,v x ,v y ) = (0.10, 1.00, 0.00, 0.99) at 0.5 < x < 1.0, 0.5 < y < 1.0, 



where Vi (i = x, y) are the i-th components of velocities and the adiabatic index is 7 = 5/3. Figure 
contours of log^p) at t = 0.4. 



shows color 




0.2 0.4 0.6 



Fig. A2. — Two-dimensional shock tube problem at t = 0.4 with 10 3 X 10 3 square uniform meshes. Color contours of the logarithm of 
the rest mass density are plotted. 
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Double Mach reflection problem 

IWoodward fc Colellal (|1984f ) introduced a double Mach reflec tion problem of a strong sh ock in the Newtonian case 
and the problem was extended to the special rel ativistic case ( Zhang & MacFadyen 2006J) . I te st the code with the 
same initial conditions as in the previous studies (jZhang fc MacFadv cn 2006; Mi zuta et al.ll2006f h The computational 
region isa0<ir<4, 0<y<l rectangle captured by 512 x 128 uniform square meshes. The density and pressure of 
the unshocked gas, the adiabatic index, and the classical Mach number (M c = Vs/c s , where Vs is a shock velocity and 
c s is a sound speed) are set to be p = 1.4, p — 0.0025, 7 = 1.4, and M c = 10 respectively. Thus, the shock velocity is 
Vs = 0.4984 and the density, pressure, and velocity of the post-shock gas are p = 8.564, p = 0.3804, and v = 0.4247. 
The shock front is initially set to cross the x axis at x — 1/6 with an angle of 60°. The initial conditions are as follows: 



(p,p,v x ,v y ) = (1.4, 0.0025, 0.00, 0.00) at x -y tan 60° < 1/6, 
(fi,p,v x ,v y ) = (8.564, 0.3804, 0.4247 sin60°, 0.4247 cos60°) at x 



y tan 60° > 1/6. 



The boundary conditions are inflow of the post-shock gas at x — 0, at < x < 1/6 on y — axis, and at < x < 
(1/6 + tan 60° + Vgt/ sin 60°) on y — 1 axis, reflect ive b oundaries at 1/6 <x<Aony — axis, and freely inflow and 
outflow conditions at the other boundaries. Figure [A3] shows thirty-level iso-surfaces of log(p) at t = 4.0. 




Fig. A3. — Double-Mach reflection problem at t = 4.0 with 512 X 128 uniform square meshes. Thirty equally-spaced contours of the 
logarithm of the rest mass density are plotted. 



Emery step problem 

A wind tunnel problem with a re flecting step, called the Emery step problem, was introduced by lEmervi Jl968) 
and repeated by subsequent studies (| Zhang fc MacFad vcn 20061; [Mizuta et al]|2006[ ). The computational domain is a 
0<ir<3,0<y<l rectangle captured by 240 x 80 uniform square meshes. The domain is filled by a gas with p = 1.4, 
p = 0.1534, and v = 0.999. The pressure is set so that M c — 3. The wind tunnel has the reflecting step locating at 
0.6 < x and y < 0.2. The following boundary conditions arc applied: the upper and lower boundaries are reflective 
and the boundaries at x = and x — 3 are open. Figure lA"4l shows thirty-level iso-surfaces of log(p) at t — 4.0. 




Fig. A4. — Emery step problem at t = 4.0 with 240 X 80 uniform square meshes. Thirty equally-spaced contours of the logarithm of the 
rest mass density are plotted. 



NUCLEOSYNTHESIS CALCULATION 

The two-dimensional special relativistic hydrodynamical calculation does not include nuclear energy releases and 
applies a constant adiabatic index 7 = 4/3. The proper internal energy is written as e pr 4 n = p/(7 — 1). Temperat ure 
T is derived with an analytical equation of state including radiation and e~-c + pair (e.g.. lFreiburghaus et al.|[T999D as 
follows: 

7 



--pr,m 



aT^ll 



T 2 



where a = 7.57 x 10 ergs cm 



4 T 9 2 + 5.3 
-3 j^-4 j g ^g ra di a tion-density constant and Tg 



(Bl) 



T/10 9 K. 
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I calculate explosive nucleosynthesis in a non-relativistic spherical explosion of a 40Af Q star with 3 x 10 52 ergs. 
The post-processing calculations are performed with the thermodynamic histories obtained by the two-dimensional 
special relativistic Eulerian hyd rody namics code and a one- dimensional Lagrangian hydrodynamics code used in 
lUmeda fc Nomotol (|2002t [2005) and iTominaga et al.l (|2007bD . The one-dimensional Lagrangian hydrodynamical 
calculation includ es nuclear energy releases from the a- network and the equation of state includes the gas, radia- 
tion, e _ -e + pair (Sugimoto & Nomoto 1975), Coulomb interactions between ions and electrons, and phase transition 
(|Nomoto et al. 19 82; N omoto fc Hasmmotol ll988) . The abundance distributions derived by both calculations are con- 
sistent (Figs. IBla b). This consistency justifies the results of the two-dimensional special relativistic hydrodynamics 
and nucleosynthesis calculations even with the simple assumptions. 
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Fig. Bl. — Abundance distributions after the non-relativistic spherical explosion of a 4QMq star with 3 X 10 ergs. The post-processing 
calculations are performed with the thermodynamic histories obtained by (a) the one-dimensional spherical hydrodynamics code and (b) 
the two-dimensional special relativistic hydrodynamical code. 



JET INJECTION 

The jet characterized by the five parameters, Edcp, E^ op , 9j ot , Tj Ct , and /th, is injected from the inner boundary at 
Ro- The jet injection is performed by putting the velocities (tv.jet, u e,jct), pressure (pjct), and density (pj e t) of the jet 
at the inner boundary as the boundary condition, iVjet, vsjct, Pjct, and pj C t are derived from the five parameters. 

The superficial area of the jet injection is Aj Ct — AttRq(1 — cos#j t) and thus the energy deposition rate is -Edcp = 
Ajet/ejcti where / e jet is an energy density flux of the jet. / e jet is written with the jet properties as / e jet == S, 



Dj ct v r ,ict, where SVjet 



PjethjetljetVrjeti ^jct — PietTiet, and h 



J r,jet 

Pj et i j et , ana ftjet is the specific enthalpy of the jet. Although the 



energy density flux at the inner boundary (/ e ) depends on the properties of the jet and the infalling matter and changes 
with time, f e is equivalent to / e j e t when the jet is injected freely. 

On the other hand, the energy density per unit volume is ej e t = Tj e t = Pjct^jctr, 2 0t — Pj e t — PjetTjct- The energy 
density consists of the thermal energy density (ethjet) and the kinetic energy density (ekinjct)- Applying /ij et 



1 + e pr ,in,jet/Pjet +Pjet/Pjet and e 



pr,in,jet 



Lorentz factor rj ct , and adiabatic index 7 of the jet as follows: 



: Pjet/il - 1): ethjet and e k i n jet are written by the density p jot , pressure p jct , 



ekinjet — (1 — /th)Sj. 
Gthjet — fth^jet = 



ft«t 



ftct (r jct 
7 



Tjet) , 
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-Ljct i 



Since 



jet = /ejet /^V jet — Pjeti the velocities, pressure, and density of the jet are written as follows: 



'■■ =^i-i/r 2 ct , 
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